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ABSTRACT 

The  behavior  of  a  fully-submerged  hydrofoil  has  been 
studied  by  using  Root  Locus  technique.   A  multivariable  con- 
cept was  used  in  designing  the  compensators  of  the  system. 
A  digital  computer  simulation  program  using  Digital  Simula- 
tion Language  or  DSL  is  produced  to  study  the  performance  on 
the  digital  computer  IBM  system  360/67.   The  Highpoint 
PC(H)-1  is  used  as  the  model  of  this  simulation. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 14 

A.  FOIL  TYPES 14 

B.  FOIL  ARRANGEMENT 16 

II.  .HYDROFOIL  EQUATIONS  STUDY  21 

A.  COORDINATE  SYSTEMS  21 

B.  EQUATION  OF  MOTION 23 

C.  CALCULATION  OF  VELOCITY  COMPONENTS  27 

D.  ANGLE  OF  ATTACK  AND  ANGLE  OF  SIDE  SLIP  30 

E.  EXPANSION  OF  FORCE  AND  MOMENT  TERMS  32 

F.  SPECIFICATION  OF  UNITS 34 

III.  HYDROFOIL  SIMULATION  35 

A.  SIMULATION  OBJECTIVE  35 

B.  COMPUTER  SIMULATION  36 

C.  AUTOMATIC  CONTROL  38 

1.  Vertical  Gyros 38 

2.  Rate  Gyros 39 

3.  Accelerometer 39 

4.  Height  Sensor 40 

IV.  TRANSFER  FUNCTION  AND  ROOT  LOCUS  DESIGN  41 

A.  PITCH-HEAVE  MOTION  TRANSFER  FUNCTION  41 

B.  ROOT  LOCUS  DESIGN  OF  FORWARD  FOIL  AUTOMATIC 
CONTROL 49 

C.  ROOT  LOCUS  DESIGN  OF  AFT  FOIL  AUTOMATIC 

CONTROL 54 


V.  MULTIVARIABLE  SYSTEM  60 

A.  BASIC  CONCEPT 60 

1.  The  Characteristic  Equation  61 

2.  Connection  Between  Single-Loop  and 
Multivariable  62 

B.  ROOT  LOCUS  DESIGN  FOR  PITCH-HEAVE  MODE 63 

C.  ROOT  LOCUS  DESIGN  FOR  ROLL  MODE 76 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY-  82 

APPENDIX  A   SUMMARY  OF  COMPLETE  EQUATIONS  84 

APPENDIX  B   PLOTS  OF  HYDRODYNAMIC  COEFFICIENTS  89 

COMPUTER  PROGRAMS  93 

LIST  OF  REFERENCES 101 

INITIAL  DISTRIBUTION  LIST  102 


6 


LIST  OF  DRAWINGS 

1  TYPE  OF  HYDROFOIL  CRAFT 14 

2  FOILS  CONFIGURATION  17 

3-A   FOIL  IN  CANARD  ARRANGEMENT 18 

3-B   FOIL  IN  AIRPLANE  ARRANGEMENT 18 

4-A   DRAWING  OF  FOIL  AND  STRUT  CONFIGURATION 19 

4-B   DRAWING  OF  FOIL  AND  FLAP  CONFIGURATION 20 

5  ORIENTATION  OF  THE  BODY  AXES  WITH  RESPECT  TO  THE 
CRAFT  AND  DIRECTIONS  OF  POSITIVE  VELOCITIES  21 

6  DEFINITION  OF  POSITIVE  DIRECTIONS  OF  EULER  ANGLE  —  26 

7  DEFINITION  OF  SYMBOLOGY  28 

8-A   ANGLE  OF  ATTACK 31 

8-B   ANGLE  OF  SIDESLIP 31 

9  BLOCK  DIAGRAM  OF  EQUATION  RELATIONSHIPS  37 

10  ROOT  LOCUS  PLOT  OF  G     50 

Pi  i 

11  ROOT  LOCUS  PLOT  OF  G     51 

P21 

12  ROOT  LOCUS  PLOT  OF  G     52 

P12 

13  ROOT  LOCUS  PLOT  OF  G     53 

P22 

14  BLOCK  DIAGRAM  SHOWING  FORWARD  FOIL  AUTOMATIC 

CONTROL 54 

15  STEP  RESPONSE  OF  THE  CRAFT  (FORWARD  FOIL  ONLY)  55 

16  BLOCK  DIAGRAM  SHOWING  AFT.  FOIL  AUTOMATIC  CONTROL  -  56 

17  STEP  RESPONSE  OF  THE  CRAFT  (AFT  FOIL  ONLY, 

INPUT  =  2.0)  57 

18  STEP  RESPONSE  OF  THE  CRAFT  (AFT  FOIL  ONLY, 

INPUT  =  0.5)  59 


7 


19  BLOCK  DIAGRAM  SHOWING  A  MULTI VARIABLE  SYSTEM 

WITH  DIAGONAL  MATRIX  CASCADE  COMPENSATOR  61 

20  ROOT  LOCUS  PLOT  OF  1+G    -G     65 

Pi  1   Cn 

21  ROOT  LOCUS  PLOT  OF  l+(det-G  /G    )-G    66 

P   P22    en 

22  ROOT  LOCUS  PLOT  OF  1+G   «G     70 

eq   c22 

23  BLOCK  DIAGRAM  SHOWING  THE  COMPENSATED  SYSTEM  71 

24  STEP  RESPONSE  OF  HEIGHT  (2  FT.  UP) 72 

25  STEP  RESPONSE  OF  HEIGHT  (2  FT.  DOWN) 73 

26  STEP  RESPONSE  OF  PITCH  ANGLE  (BOW  UP) 74 

27  STEP  RESPONSE  OF  PITCH  ANGLE  (STERN  UP) 75 

28  ROOT  LOCUS  PLOT  OF  ROLLING  TRANSFER  FUNCTION  78 

29  STEP  RESPONSE  OF  ROLL  ANGLE  (PORT  UP) 79 

30  STEP  RESPONSE  OF  ROLL  ANGLE  (STARBOARD  UP)  SO 

31  BLOCK  DIAGRAM  OF  CONTROL  SYSTEM  81 


8 


TABLE  OF  SYMBOLS 


"S 
CG 

CDF 
CDS 


'S 


DiF 


g, 


gr 


FEET2 
FEET2 

Dimension less 
Dimension less 
Dimensionless 
Dimension less 
Pounds 


F 
DiS 

Pounds 

F 
LiF 

Pounds 

F 
SiS 

Pounds 

F 
X 

Pounds 

F 

XiF 

Pounds 

F 

Pounds 

F 
YiF 

Pounds 

F 
YiS 

Pounds 

F 
Z 

Pounds 

F 
ZiF 

Pounds 

F 
ZiS 

Pounds 

gy 

Feet/Second 

Feet/Second: 


Feet/Second: 


Area  of  foil 

Area  of  strut 

Center  of  gravity 

Drag  coefficient  of  a  foil 

Drag  coefficient  of  a  strut 

Lift  coefficient 

Side  force  coefficient 

Drag  force  on  a  foil  (water  axes) 

Drag  force  on  a  strut  (water  axes) 

Lift  force  on  foil  (water  axes) 

Side  force  on  a  strut  (water  axes) 

Force  in  direction  of  body  X-axis 

Force  on  a  foil  in  body  X-axis 

Force  in  body  Y-axis 

Force  on  a  foil  in  body  Y-axis 

Force  on  a  strut  in  body  Y-axis 

Force  in  direction  of  body  Z-axis 

Force  on  a  foil  in  body  Z-direction 

Force  on  a  strut  in  body  Z-direction 

Component  of  gravity  in  body 

X-direction 

Component  of  gravity  in  body 

Y-direction 

Component  of  gravity  in  body 

Z-direction 


9 


XX 


LYY 


ZZ 


JXHS 


XiF 


^ 


iF 


LZiF 


JZiS 


M 


Feet 


Slug-Foot : 


Slug- Foot ; 


Slug-Foot : 


Foot-Pound 


Feet 


reet 


Feet 


Feet 


Feet 


Foot-Pounds 


Height  of  Height  Sensor  above 

instantaneous  water  surface 

Moment  of  inertia  about  body 

X-axis 

Moment  of  inertia  about  body 

Y-axis 

Moment  of  inertia  about  body 

Z-axis 

Moment  acting  to  produce  roll  about 

body  X-axis 

Distance  from  CG.  along  X-axis  to 

a  point  of  application  of  foil 

force 


Distance  from  CG.  along  X- 


p  *v  "i  o   +■  r~s 


a  point  of  foil  force 

Distance  from  CG.  along  body  Y-axis 

to  a  point  of  application  of  foil 

force 

Distance  from  CG.  along  body  Z-axis 

to  a  point  of  application  of  foil 

force 

Distance  from  CG.  along  body  Z-axis 

to  a  point  of  application  of  strut 

force 

Moment  acting  to  provide  pitch  about 

body  X-axis 


10 


m      Slugs 

N      Foot-Pounds 


P 

Q 

R 
S 


EiF 


EiS 


SHP 
T 

u 

UE 

U. 

l 


V, 


w 


WE 

w. 

1 


Radians /Second 

Radians /Second 

Radians/Second 

Feet 

Feet 


Horsepower 
Pounds 

Feet/Second 
Feet/Second 
Feet/Second 

Feet/Second 

Feet/Second 
Feet/Second 

Feet/Second 

Feet/Second 
Feet/Second 


Mass  of  the  craft 

Moment  acting  to  produce  yaw  about 

body  Z-axis 

Roll  rate  about  body  X-axis 

Pitch  rate  about  body  Y-axis 

Yaw  rate  about  body  Z-axis 

Submergence  of  a  foil  in  earth  axes 

Submergence  in  earth  axes  of  the 

assumed  point  of  application  of 

strut  force 

Shaft  horsepower 

Magnitude  of  thrust  in  direction  of 

body  X-axis 

Velocity  in  direction  of  body  X-axis 

Velocity  along  earth  X-axis 

Foil  velocity  in  direction  of  body 

X-axis 

Velocity  in  direction  of  body 

Y-axis 

Velocity  along  earth  Y-axis 

Foil  velocity  in  direction  of  body 

Y-axis 

Velocity  in  direction  of  body 

Z-axis 

Velocity  along  earth  Z-axis 

Foil  velocity  in  direction  of  body 

Z-axis 


11 


E 


E 


a. 


Feet 


Feet 


Radians 


Radians 


n 

Dimension less 

p 

Slugs/Foot3 

* 

Radians 

* 

Radians 

e 

Radians 

Distance  along  earth  X-axis  from 

reference  origin 

Distance  along  earth  Y-axis  from 

reference  origin 

Angle  of  attack  of  a  particular 

foil,  designated  by  subscript 

Angle  of  side  slip  of  particular 

strut,  designated  by  subscript 

Propulsion  efficiency 

Water  density 

Roll  angle 

Yaw  angle 

Pitch  angle 


Subscripts 

C      Refers  to  Center  foil  or  strut 

P      Refers  to  Port  foil  or  strut 

S      Refers  to  Starboard  foil  or  strut 

M      Refers  to  Mid  foil  section 

F      Refers  to  foil 

S      Refers  to  strut 

i      When  "i"  appears  in  the  subscript  to  a  symbol,  it 
indicates  that  the  symbol  is  to  be  repeated  with  i 
successfully  replaced  by  C,  P,  S,  and  M,  as  indicated. 


12 


ACKNOWLEDGEMENT 

I  would  like  to  express  my  deepest  appreciation  for  Dr. 
George  J.  Thaler's  considered  advice  and  help.  For  without 
them,  this  research  paper  would  not  have  been  possible. 


13 


INTRODUCTION 


A.   FOIL  TYPES 

The  main  idea  of  the  hydrofoil  craft  is  to  reduce  the 
resistance  of  the  craft  by  reducing  the  area  in  contact  with 
the  water  and  thus  increase  the  speed  for  the  same  power. 
One  way  to  accomplish  this  is  fitting  foils  below  the  hull 
of  a  craft,  so  that  the  craft  is  lifted  clear  of  the  water 
at  high  speed.   There  are  two  principal  foil  arrangements  in 
use  today,  classed  as  surface  piercing  foil  and  fully-sub- 
merged foil,  as  shown  in  Figure  1. 


Surface-Piercing  Foil 

FIGURE  1 
TYPE  OF  HYDROFOIL  CRAFT 


Fully-Submerged  Foil 


The  surface-piercing  foil  is  supported  by  vertical 
struts  and  the  lifting  surface  extends  in  a  V-shape  from  the 
water  surface  down  through  the  water  to  the  point  of  conver- 
gence of  the  V.   The  surface-piercing  foil  provides  a 
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self-stabilization  for  the  craft  in  pitch,  height  and  roll. 
The  self-stabilizing  aspect  of  the  surface  piercing  foil  is 
easy  to  see.   Suppose  the  craft  pitches  up  slightly,  then 
less  area  of  the  forward  foil  and  more  area  of  the  after  foil 
is  submerged.   This  results  in  less  lift  from  the  forward 
foil  and  more  lift  from  the  after  foil  so  the  craft  will 
settle  back  again,  the  same  stabilization  occurs  in  roll. 
For  the  height  as  the  craft  speed  increases  the  whole  craft 
rises  up  until  the  amount  of  wetted  surface  of  the  foils  is 
just  sufficient  to  create  a  lift  equal  to  the  weight  of  the 
craft. 

The  fully-submerged  foil  cannot  adjust  itself  to  chang- 
ing lift  demands  by  a  change  in  area,  but  only  by  a  change 
in  incidence,  either  of  the  whole  foil  or  of  attached  flaps. 
This  involves  the  use  of  some  sensing  arrangement  to  control 
the  incidence  by  monitoring  the  height  of  the  craft  above 
water  surface.   The  fully-submerged  foil  is  a  more  sophisti- 
cated design,  and  the  control  mechanism  is  likely  to  be  com- 
plicated.  But  the  fully-submerged  foil  has  a  very  distinct 
advantage  when  operating  in  heavy  seas.   The  surface-piercing 
foil  will  respond  to  even  relatively  small  changes  in  the 
surface  of  the  sea  and  give  a  very  rough  ride.   The  fully- 
submerged  foil,  on  the  other  hand,  is  much  less  affected  by 
the  change  in  the  surface  of  the  sea.   This  paper  will  deal 
only  with  the  fully-submerged  foil  craft. 
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B.   FOIL  ARRANGEMENT 

Most  of  the  fully-submerged  foil  craft  have  two  foils, 
one  at  the  forward  and  the  other  at  the  aft.   Both  foils  are 
supported  by  the  vertical  struts.   Figure  2  shows  the  con- 
figuration of  main  foil  and  secondary  foil;  Figures  3A  and 
3B  show  foils  in  canard  arrangement  and  in  airplane  arrange- 
ment respectively. 

The  simulation  in  this  paper  is  based  on  the  canard  ar- 
rangement, which  has  a  secondary  foil  as  a  forward  foil  and 
a  main  foil  as  the  aft  foil.   The  HIGH  POINT  PC(H)-1  was  the 
model  for  this  simulation.   Drawing  of  the  craft  and  the 
foils  are  shown  in  Figures  4A  and  4B. 
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FIGURE  2 
FOILS  CONFIGURATION 
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FIGURE  3A 
FOIL  IN  CANARD  ARRANGEMENT 


FIGURE  3B 
FOIL  IN  AIRPLANE  ARRANGEMENT 
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II.   HYDROFOIL  EQUATIONS  STUDY 

A.   COORDINATE  SYSTEMS 

Three  coordinate  systems  will  be  used  in  this  study  of 
hydrofoil  craft  dynamics.   They  are:  the  Body  Axes,  the 
Water  Axes  and  the  Earth  Axes.   Each  is  in  a  right  hand 
orthogonal  system. 

The  body  axes  coordinate  system  (x,y,z)  shown  in  Figure 
5  has  its  origin  at  the  center  of  gravity  of  the  craft  and 
is  fixed  relative  to  the  craft.   This  system  is  initially 
coincident  with  the  origin  of  the  earth  axes  system,  however, 
it  moves  with  the  craft  as  time  progresses.   X  is  positive 
forward,  Y  is  positive  to  starboard  and  Z  is  positive  down- 
ward. 


FIGURE  5 
ORIENTATION  OF  THE  BODY  AXES  WITH  RESPECT  TO  THE 
CRAFT  AND  DIRECTIONS  OF  POSITIVE  VELOCITIES 
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The  earth  axes  coordinate  system  (X^.Y-.Z^)  is  fixed 
relative  to  the  earth  surface.   The  origin  may  be  chosen  at 
any  convenient  location,  provided  the  body  system  origin  is 
initially  at  the  same  point.   The  X„  axis  lies  in  the  hori- 
zontal  plane;  it  is  initially  coincident  with,  and  positive 
in  the  same  direction  as  the  body  X  axis.   The  Yp  axis  lies 
in  the  horizontal  plane  and  is  positive  to  starboard  when 
the  observer  is  facing  the  positive  XF  direction.   The  Z_ 
axis  is  normal  to  the  horizontal  and  is  positive  downward. 
Water  surface  motion  and  craft  motion  relative  to  calm  water 
are  described  in  this  coordinate  system. 

The  water  axes  coordinate  system  is  aligned  with  the 
relative  velocity  vector  and  resolves  into  "Lift,"  "Drag," 
and  Side -Force"  direction. 


DRAG 

It  will  be  specified  that  the  order  of  rotation  from 
water  axes  to  body  axes  is: 

1)  The  (D,S)  axes  rotate  about  the  L-axis  through  the 
angle  3  to  the  position  (D_.,S  )  where  S1  is  coincident  with 
the  body  Y-axis. 

2)  The  (D  ,L)  axes  rotate  about  the  S.-axis,  or  more 
pertinently  about  the  Y-axis,  through  the  angle  a  to  the 
position  (D   L  )  which  corresponds  to  the  (-X.-Z)  axes. 
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B.   EQUATION  OF  MOTION 

The  development  of  the  six  equations  of  motion  for  a 
hydrofoil  follows  the  same  procedure  as  those  for  a  displace- 
ment type  hull.   This  development  is  clearly  outlined  in 
Ref.  2  and  8  and  therefore  will  not  be  done  here. 

The  equations  of  motion  are: 
Translation 

u  =  -Fv  -  QW  +  RV  (2-1) 

in  a 

v  =  -Fv  +  PW  -  RU  (2-2) 

m  Y  v    ' 

™  =  mFZ   +   QU   "   PV  (2_3) 

Rotation 

P   =  y^-[L-QR(Izz-Iyy)+(R+QP)Ixz]  (2-4) 

Q   =  t^-[M-RP(Ixx-Izz)-(P2-R2)Ixz]  (2-5) 

R   =    j^-[N-PQ(IYY-Ixx)-(QR-P)Ixz]  (2-6) 

ZiZt 

The  terms  RV,  PV ,  and  PW  are  very  small  and  may  be  neg- 
lected.  Similarly,  the  product  of  inertia  quantity,  I,r7>  is 
less  than  10%  of  I„x  and  less  than  4%  of  I77-   Consequently, 
when  considered  in  conjunction  with  the  magnitudes  of  (R+PQ), 

XZ 


(P2-R2 )  and  (QR-P),  the  terms  containing  IY7  may  be  neglected 


for  this  craft. 

Equations  (2-1)  through  (2-6)  are  linear  and  rotational 
accelerations  in  body  axes.   The  velocity  terms  in  body  axe. 
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can  be  obtained  directly  by  integrating  each  equation  one 

time. 

U  =   /  Udt  (2-7) 

V  =   /  Vdt  (2-8) 

W  =   /  Wdt  (2-9) 

P  =   /   Pdt  (2-10) 

Q  =   /  Qdt  (2-11) 

R  =    /   Rdt  (2-12) 

Linear   displacements   and   angular  positions  with  respect 
to  body   axes   could  be  obtained  by   an   additional   integration, 
but   such  quantities   are  of    limited  usefulness.      The   dis- 
placements  and   angle  usually   of   interest   are  those  with 
respect    to   earth   axes.       In   order   to   calculate  the   craft's 
position   in   earth   axes   coordinates,    the   craft   velocity   vec- 
tor must   be   transformed   from  body   axes   coordinates    into 
earth  axes   coordinates   as   follow: 

U,-,   =   Ucos6cosij>+V[cosijjsin0sinij;-Sinifjcos<f>] 

+W[cosijjsin0cos<j)+sin<j>sincj)]  (2-13) 

Vp   =   Ucos0sinip+V[cosijjcos<j)+Sinijjsin0sin<f>] 

+W[sinij;sin0cos(f>-cos^sinci)]  (2-14) 

W£  =  -Usin8+Vcosesin<t>+Wcos9cos(J>  (2-15) 

Only  the  craft's  velocity  in  the  downward,  ZF,  direction 
is  required  to  calculate  its  height  above  the  water  surface. 
That    is,    the   vertical    position    in    earth   axes    follows    from 

ZE  =  /  WEdt  (2-16) 
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The  craft's  position  coordinates  are  obtained  from 

X£  =  /  UEdt  (2-17) 

YE  =  /  V£dt  (2-18) 

The  Euler  angles,  <J> ,  8,  and  iJj  are  also  of  interest  and 

•     •  • 

can  be  found  by  integration  of  Euler  rates  <}> ,  0,  and  ty .  The 
rates  in  turn  must  be  derived  from  the  body  angular  rates. 
The  Euler  rates  are  not  easily  obtained  because  they  occur 
about  axes  which  are  not  orthogonal.  This  fact  can  be  ap- 
preciated by  referring  to  Figure  6  and  recalling  how  the 
angles  were  defined.  f  is  a  rotation  about  the  Z0  axis,  0 
is  a  rotation  about  the  Yi  axis,  and  f  is  a  rotation  about 
the  X  axis.   They  result  in  the  following  set  of  equations: 

$  =  p+<jjsin0  (2-19) 

0    =   Qcoscj>-Rsin<j)  (2-10) 

4>  =  (Qsin<(>+Rcos(|>)cos6  +  ($-P)sine  (2-21) 

Studies  of  this  craft  show  that  even  under  the  most 
drastic  abnormal  situations,  pitch  angles  are  not  expected 
to  reach  10°  and  under  normal  operating  conditions,  roll 
angles  will  not  exceed  3°.   Therefore,  the  following  assump- 
tions are  justified; 

sine  =  0 

(2-22) 

COS0  =  1 

The  Euler  angles  are  obtained  by  integration  of  equations 
(2-19)  through  (2-21). 
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FIGURE  6 
DEFINITION   OF   POSITIVE   DIRECTIONS    OF    EULER   ANGLE 
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C.   CALCULATION  OF  VELOCITY  COMPONENTS 

The  total  velocity  components  of  each  foil  in  body  axes 
can  be  expressed  in  terms  of  craft  linear  and  angular  veloc- 
ities by  considering  the  craft  geometry  as  represented  in 

Figure  7  and  Table  I. 

The  resultant  equations  for  the  example  shown  are: 


Center  Foil 


^C  -  W  "  LXCF* Q 


Port  Foil 


VS  -  V  -  LZSF  P  +  LXSF  R 


Wg  =  W  +  L^P  -  ^XSF 
Mid-Foil  Segment 

UM  =  U  +  LZMF  Q 


\   =  W  "  LXMF  Q 


(2-23) 


UC  =  U  +  LZCF*Q 

VC  =  V  -  LZCF'P  +  LXCF-R  (2-24) 


(2-25) 


-   Up  =  U  +  LZPF*Q"  W*  (2_26) 

VP  =  V-  LzpF  P  +  Lxp/  R  (2~27> 

Wp-  W  +  LYpF'P  -  LxpF'Q  <2-28> 
Starboard  Foil 

«S  -   »  +  LZSF  Q  "  LVSF-  R  ^'^ 


(2-30) 


p  -  lwQ  (2"31) 


(2-32) 


VM  =  V  -  LMF- P  +  LXMp- B  C2-33) 


(2-34) 
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FIGURE  7 
DEFINITION  OF  SYMBOLOGY 

o  =  Point  of  application  of  foil  forces 
o  =  Point  of  application  of  strut  forces 

LXCF  =  LXCS 

LXPF  "  LXPS  "  LXSF  =  LXSS  =  LXMF 

^sf  =  Srss  :  lypf  =  -Sps 

LZPF  =  LZSF  =  LZMF 

Lv^o  and  L  cc  are  variables  and  in  general  will  not  be 

irb  Zoo         , 

equal . 
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These  foil  velocities  are  in  body  axis  terms  and  can  be 
used  directly  to  calculate  the  angles  of  attack  and  angles 
of  side  slip. 

D.   ANGLE  OF  ATTACK  AND  ANGLE  OF  SIDE  SLIP 

The  angle  of  attack  of  the  foil  is  shown  in  Figure  8A 
and  defined  as  the  angle  whose  tangent  is  the  total  relative 
velocity  between  the  water  and  the  foil  in  the  body  Z  direc- 
tion, divided  by  the  total  relative  velocity  between  the 
water  and  foil  in  the  body  X  direction 

W 
a  =   arc  tan  =—  (2-35) 

r 

If  the  foil  is  not  aligned  with  the  body  X-axis,  then 
the  foil-  has  a  fixed  angle  of  attack  and  the  total  angle  of 
attack  becomes 

aTOTAL  =  aFIXED  +  a  (2_36) 

The  side  slip  angle  is  shown  in  Figure  SB  and  can  be 
represented  by 

Vri 

3.  =  arc  tan  -TZZZZZZZZZZ^  (2-37) 

1         /u2.+v2.+w2. 

r  1   r  1   r  1 

During  actual  craft  operations,  the  angles  will  be  less 
than  10°,  and  will  permit  use  of  the  small  angle  approxima- 
tion.  Also  the  V  and  W  velocities  can  be  considered  very 
small  in  relation  to  the  U  velocity.   Applying  these  simpli- 
fications, the  equations  for  a  and  6  become: 
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E.   EXPANSION  OF  FORCE  AND  MOMENT  TERMS 

With  the  aid  of  Figure  7  and  Table  I,  the  forces  and 
moments  are  set  down.   It  is  important  that  when  a  numerical 
value  is  inserted  for  one  of  the  alphabetical  dimensions, 
the  sign  of  the  value  must  be  positive  or  negative  according 
to  positive  and  negative  directions  from  the  axis  system 
origin  at  the  C.G.   Also  note  that  horizontal  foils  make  no 
contribution  to  the  body  Y  force  and  vertical  struts  do  not 
contribute  to  the  body  Z  force.   One  of  the  assumptions  made 
upon  commencing  simulation  was  that  the  struts  contributed 
zero  lift  and  that  the  foils  contributed  zero  side  force. 
Therefore,  all  body  Y  motion  originates  at  the  struts  and 
all  body  Z  motion  originates  at  the  foils. 

The  force  and  moment  equations  become: 

FX  =  FXCF+FXPF+FXSF+FXCS+FXSS+mgx+TX  (2~38) 

FY  =  FYCS+FYPS+FYSS+mgy  ■    (2"39) 

FZ  =  FZCF+FZPF+FZSF+FZMF+mgz  (2"40) 

L   =  (FzLY)pF+(FzLY)SF-(FYLz)ps-(FyLz)ss-(FYLz)cs      (2-11) 

M   =  -  (  FZLX )  c]p-  ( FZLX )  pF-  ( FZLX  )  sp-  ( F^ )  Mp+  ( F^  )Q¥ 

+(FxLz)pF+(FxLz)SF+(FxLz)MF+(FxLz)cs+(FxLz)ps 

+  <FXLZ>SS+TXLZT  (2~'12) 

N  =  ^y4)Cs+<pylk)ps+(py4)ss  (2~43) 
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All  of  these  forces  have  a  hydrodynamic  origin  except 
for  the  thrust  and  gravity  terms.   The  thrust  is  always 
associated  with  the  body  axis  and  requires  no  further  expan- 
sion.  The  gravity  terms  can  be  simplified  using  small  angle 
approximation  for  the  Euler  angles  (see  equations  A-31,32,33) 

The  general  expression  for  hydrodynamic  forces  in  water 
axes  coordinates  is 

F.  =  qAC.  i  =  L,D,S  (2-44) 

where  L,D,S  represent  lift,  drag  and  side  force  respectively. 
The  lift,  drag,  and  side  force  coefficients  very  as  func- 
tions of  angle  of  attack,  angle  of  side  slip,  submergence, 
velocity,  flap  elevator  and  rudder  positions. 

Ideally,  mathematical  expressions  would  have  been  de- 
veloped to  correctly  depict  the  interrelation  of  all  the 
variables  which  affect  the  force  coefficients.   However, 
hydrofoil  technology  has  not  advanced  to  a  position  which 
would  yield  such  expressions.   Some  of  the  major  problems 
encountered  in  deriving  a  mathematical  expression  are: 

1)  The  occurrence  of  cavitation  and  ventilation  which 
are  completely  unpredictable. 

2)  Even  more  unpredictable  is  the  cessation  of  cavita- 
tion and  ventilation. 

3)  The  nonlinearity  of  the  hydrodynamic  coefficients. 

4)  Lack  of  sufficient  test  data  to  completely  describe 
all  of  the  above. 
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F.   SPECIFICATION  OF  UNITS 

Considerable  confusion  resulted  from  the  interpretation 
of  the  various  references  concerning  the  units  used  in  de- 
fining Euler  angles  and  angular  rotations.   This  problem  has 
been  alleviated  here  by  the  specific  definitions  in  the 
Table  of  Symbols. 
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III.   HYDROFOIL  SIMULATION 

A.   SIMULATION  OBJECTIVES 

The  object  of  the  simulation  was  to  commence  the  runs 
with  the  craft  foilborne  and  in  a  steady  state  condition. 
Once  this  condition  was  achieved,  the  step  response  of  the 
craft  would  be  studied  by  applying  step  functions  to  excite 
motion  along  and  about  each  body  axis.   That  is,  perturba- 
tions should  be  introduced  to  separately  excite: 

1)  Motion  along  the  X-axis 

2)  Motion  about  the  X-axis 

3)  ,  Motion  along  the  Y-axis 

4)  Motion  about  the  Y-axis 

5)  Motion  along  the  Z-axis 

6)  Motion  about  the  Z-axis 

After  the  step  response  was  obtained,  a  sinusoidal  sea 
would  be  inserted  to  observe  the  craft's  motion  in  a  seaway. 
The  following  assumptions  were  made  prior  to  beginning  simu- 
lation : 

1)  Craft  equations  of  motion  are  valid  only  for  the 
foilborne  mode. 

2)  Weight  of  the  craft  remains  constnat. 

3)  Hydrodynamic  coefficients  are  based  on  fully  wetted 
surfaces,  i.e.  no  cavitation  or  ventilation. 

4)  No  constraint  was  placed  on  craft's  heading. 
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5)  Foils  are  considered  flat  surfaces  vice  dihedral/ 
anhedral. 

6)  The  craft  may  be  operated  in  calm  water  or  seaway. 

B.   COMPUTER  SIMULATION 

The  computer  language  used  in  this  simulation  was  Digital 
Simulation  Language  or  DSL,  and  based  on  the  computer  pro- 
gram in  Ref.  1. 

Before  the  main  program  could  be  assembled,  subroutines 
had  to  be  written  to  find  the  hydrodynamic  coefficients  of 
foil  lift  and  drag,  and  strut  drag  and  side  slip.   Using  the 
simulation  curves  in  Ref.  3  as  a  guide,  the  curves  were 
stored  as  data  points  in  the  setup  routine.   The  subroutines 
"INTERP"  and  "INTERP1"  are  interpolation  routines  to  obtain 
the  proper  coefficient  from  the  curves.   For  INTERP,  values 
of  angle  of  attack  and  flap  angle  are  used  to  enter  the 
curves  and  obtain  a  value  of  foil  drag.   Angle  of  side  slip 
and  submergence  are  used  to  enter  INTERP1  to  obtain  strut 
drag.   The  curve  for  the  lift  coefficient  of  a  foil  in  the 
fully  wetted  region  is  straight  line  so  this  coefficient  is 
found  by  merely  solving  the  equation  for  the  straight  line. 
Curves  plotted  from  the  subroutines  are  shown  in  Appendix  B. 
Listings  of  the  subroutines  are  shown  in  the  computer  pro- 
gram section.   The  main  program  was  then  assembled  using 
Figure  9  as  a  reference  for  data  flow.   One  part  before  the 
end  of  the  main  program  was  opened  for  an  automatic  control 
block.   A  complete  listing  of  the  main  program  is  shown  in 
the  computer  program  section. 
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C.   AUTOMATIC  CONTROL 

Control  signals  originate  at  the  pilot's  controls,  at  the 
motion  sensors,  and  at  the  position  transducer.   These  are 
combined  and  processed  in  the  electronics  to  produce  the 
signal  for  the  servo-valve.   The  pilot  will  be  able  to  move 
the  helm,  a  lever,  or  a  knob  to  introduce  steering,  altitude 
and  attitude  commands,  the  motion  sensors  sense  errors  be- 
tween the  commanded  and  actual  craft  motions  and  produce 
electrical  signals  proportional  to  these  errors.   The  posi- 
tion transducer  provides  the  feedback  signal  so  that  the 
control  surface  stope  moving  when  it  has  reached  the  position 
corresponding  to  the  processed  error  signals. 

It  is  necessary  to  consider  how  the  control  system  com- 
ponents affect  the  equations  cf  motion  through  the  control 
surface  deflections.   In  general,  a  signal  originating  at 
the  pilot  or  at  the  motion  sensors  will  be  modified  by  the 
individual  dynamics  of  the  components  in  the  signal  path  be- 
fore the  control  surface  motion  actually  occurs.   Motion 
sensors  to  be  discussed  are  vertical  gyros,  rate  gyros,  ac- 
celeromaters  and  height  sensors. 

1.   Vertical  Gyros 

A  vertical  gyro  of  the  type  of  interest  here  will 
produce  one  signal  directly  proportional  to  the  pitch  angle 
of  the  craft  and  one  signal  directly  proportional  to  the 
roll  angle  of  the  craft. 
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2.  Rate  Gyros 

A  rate  gyro  will  sense  an  angular  rate  about  one  of 
the  three  body  axes.   (Three-gyro  "rate  packages"  are  avail- 
able if  required  for  a  particular  application.)   The  dynam- 
ics of  the  gyro  are  such  as  to  produce  an  output  signal 
which  is  related  to  the  actual  angular  rate  by  a  second  order 
differential  equation.   Using  the  roll  rate  P,  as  an  example, 
the  equation  is 

P  "  ^?(0  +  ««»  I  +  ■»**>  (3-D 

n 

or  in  transfer  function  form  for  all  rotational  rates: 

P  _  Q  _   R  _  ^ ,«   9, 

P        Q   ~   R   "    S2+2C03    S+o)2  ^   A) 

n    n 

where 

P,Q,R  =  Outputs  of  rate  gyros  (RAD/SEC) 

P,Q,R  =  Craft  rotational  rates  (RAD/SEC) 

a)     =  Natural  frequency  of  a  particular  gyro 

(RAD /SEC) 
C     =  Damping  ratio  of  a  particular  gyro 

(Dimension less) 
S     =  Laplace  transform  variable 

3.  Accelerometer 

Outputs  of  accelerometers  are  also  assumed  to  relate 
to  the  craft  motions  by  second  order  differential  equations. 
In  general  form,  the  equation  is: 

1    d2a.        da.. 

a,  =  — ( }  +   2cu)  -r-i  +  w2a.)  (3-3) 

i    a)2   dt2       n  dt      n  l 
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or: 

h  < 

77  ~   S2+2Coj  S+o)2  (3_4) 

i         n    n 

where  i  will  become  a  symbol  corresponding  to  the  location 
of  a  particular  accelerometer . 

To  identify  the  accelerations  sensed  by  the  instru- 
ments, it  will  be  assumed  that  there  is  one  accelerometer  at 
the  C.G.  mounted  such  as  to  measure  lateral  acceleration, 
and  one  accelerometer  above  each  foil  mounted  in  the  X-Y 
plane  so  as  to  measure  accelerations  parallel  to  the  Z-axis. 

The  lateral  acceleration  at  the  C.G.  is: 

F 

a  =  -2  (3-5) 

y    m 

Taking  angular  acceleration  into  account,  the  accel- 
erations in  the  Z-direction  at  points  above  a  center,  port, 
and  starboard  foil  are: 

Azc  =  IT  "  lxcf'Q  <3-6> 

AZP  =  it   +   LXPF*Q  "  LYPF^  (3"7) 

AZS  =  -F  +  LXSF'Q  +  LYSF^  '         ■    <3~8) 

4.   Height  Sensor 

Ultrasonic  and  sonar  types  are  either  in  use  or  ex- 
pected to  be  used.   For  this  study  it  will  be  assumed  that 
the  dynamic  characteristics  of  the  sensor,  whatever  the 
type,  will  be  such  as  to  give  a  direct  indication  of  the  in- 
stantaneous height  of  the  sensor  above  the  water. 
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IV.   TRANSFER  FUNCTION  AND  ROOT  LOCUS  DESIGN 

A.   PITCH- HEAVE  MOTION  TRANSFER  FUNCTION 

The  best  way  to  study  the  stability  of  the  craft  was  to 
find  the  transfer  functions  of  the  craft  and  use  either  Root 
Locus  plots  or  Bode  plots  to  be  an  aid  in  a  stability  study 
and  also  in  designing  the  compensators  of  the  system. 

The  transfer  functions  of  the  craft  can  be  derived  di- 
rectly from  the  equations  of  motion.   In  the  pitch-heave 
mode,  equations  (2-3)  and  (2-5)  were  the  equations  that  de- 
scribed the  motion  in  the  Z-axis  and  motion  about  the  Y-axis, 
After  neglecting  the  small  terms,  equations  (2-3)  and  (2-5) 
were  as  follow 

W  =  ^*FZ+Q-U  (3-9) 

Q  =  y^-  -M  (3-10) 

YY 

From  equations  (2-35)  and  (2-36)  the  angle  of  attack  of 
each  foil  can  be  written  as  follows: 

WRC 
ttC  =  "IT  +  ALF  (3-H) 

w 

ag  =  -[p  +  ALg  (3-12) 

W 

aP  =  "IT  +  ALP  (3-13) 

w 

aM  =  "IT  +  alm  (3~14) 
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WRP  =  W  +  ^PF*1'  -  LXPF*Q  (3"15) 

WRS  =  W  +  LYSF* P  "  hiSF   Q  (3-16> 

WRM=  W+  LXMF*Q  (3-1?) 

wrc  =  w-lxcf'q  (3"18) 

Assume  there  is  no  significant  coupling  from  the  roll, 
yaw,  sway  or  surge  motions.   Equations  (3-15),  (3-16),  (3-17) 
and  (3-18)  will  be 

WRP  =  W  "  LXPF*Q 

WRS  =  W  "  LXSF*  9 

WRM  =  W  "  LXMF'Q 

WRC  =  W  "  LXCF*Q 
And  equations  (3-11),  (3-12),  (3-13)  and  (3-14)  will  be 

_  W    LXCF      flT 
aC  "  U  "  ~U Q  +  ALF 

W    LXSF  n  .  AT 

as  ~  u  "  "u — Q  +  ALS 

W    LXMF  n    AT 
aM  -  U  "  "U Q  +  ALM 

W    LXPF  n  ,  AT 

ap  "  u  "  ~u — Q  +  alp 

^SF  =  LXPF  =  LXMF 
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Let  ALp  =  Al ,    ALg  =  ALp  =  ALM  =  A2 


'XCF-Q  =  Qi 


u 

^SF 
U 


•Q  =  Q2 


W      ' 

—  =  w 
u 


ac  =  W  -Q1+A1 


ag  =  W  -Q2+A2  =  ap  =  «M 


From  Equation  (2-40)  Y^   =  F^+P^+F^+F^+P^ 

FZCF  =  -FDCz'aC-FLCZ  (3"19) 

FZPF  =  -FDPZ*ap-FLPZ  (3"20) 

FZSF  =  -FDSZ*aS-FLSZ  '                  (3"21) 

FZMF  =  "FDMF'aM~FLMZ  (3-22) 

Equation  (2-44)  shows  the  general  expression  for  hydro- 
dynamic  forces,  and  from  Ref .  3  the  linear  approximation  of 
hydrodynamic  coefficients  are: 

CL„„  =  5.15  a„  +  0.036  5„  +  0.13 
CF         C  F 

CL„„  =  4.10  a^  +  0.018  6n 
PF  P  P 

CLgF  =  4.10  ag  +  0.018  6g 

CL.„  =  5.15  ex.,  +  0.036  6. 

MF  M  A 
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CDF  =  °-012  +  0.0015  6F 


C^A  =  0.02  +  0.0028  6. 
DA  A 


CDP  =  °'02  +  °-0028  6p 


CDS  =  °'02  +  °-0028  6s 


Let 


0.5-p-U2-ACF  =  Di 


'  -\pF   =  0.5-P-  U2  „SF 


0.5'P'U2-A„„  =  0.5'P"  U2-Ac„  =  D2 


0.5-p-U2-AMF  =  D3 
After  dropping  the  nonlinear  terms 
FZCF  =  -(5-162Di )w'+(5.162D1 )Q1-D1(0.0015A1+0.036)6F 

-DL(5.162A1+0.13) 
FZPF  =  -(4.12D2)W'+(4.12D2)Q2-4.12D2A2 
FZSF  =  -(4.12D2)W'+(4.12D2)Q2-4.12D2A2 

FZMF  =  -(5-17D3)w'+(5.17D3)Q2-D3(0.0028A2+0.036)8A-5.17D3A2 
Fz    =  -(5.162DJ+8.24D2+5.17D3 )^  +(214. 69DX -139. 92D2-87. 79D3 )§ 

-Dj (0.0015A1+0.036)6F-D3 (0. 0028A2 +0. 036)6 A 

-[Dx (5.162A1+0.13)+8.24D2A2+5.17D3A2 ] 

"  _       (5.162D1+8.24D?+5.17D,)    a     ,     r ( 214 . 69Dj-139. 92D9-87 . 79D , ) 

L  mU  Z   +    [  mU 

xinj±   FG    .         D,f0.0015A! +0.036).         D,  (0.  0028AC)+0.  036  )  , 

+U16+  —   6    -   — L-^ -6-,-,-  — " i ^o, 

J  m  m  F  m  A 

[Di (5.162A,+0.13)+8.24D,A?+5.17D,A? }  (3-23) 

m 
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M     =  _M    -M    +F    L    +F    L    +F    L    +F    L    +F    L 
"         FWD   AFT   XCF  ZCF   XSF^ZSF   XMF  ZMF   XCS  ZCS   XPS  ZPS 

+FXSSLZSS+TXLZT 

M     =  F    L 
FWD     ZCF  XCF 

MAFT  ~  (FZPF+FZSF+FZMF)LXMF 

M     =  [(214.69D1)w'-(214.69D1)Q1+(0.212D1A1+1.497D1)6 
+  (214-69D1A1+5.41D1 ) ] +[ -(139. 92D2)W +( 139 . 92D2 )Q2 
-139.92D2A2-(87.79D3 )W* +(87. 79D3 )Q2-(0. 048D3 A2+0. 611D3 )6. 

-87.79D3A2]+[(281.1D2A2)w'-(281.1D2A2)Q2+(140.55D2Ai 

-0.686D2)]+[ (159.65DiAi+2.02Di )W* -(159 . 65Di Aj+2 . 02Di )Q! 

+  (0.558DiAi-0.023D! )6p+(79 . 83D: A? +2 .  C^AiDj -0 . 186) ] 

+  E(176.5D3A2)Y.T'-(176.5D3A2)Q2  +  (0.62D3A2-0.048D3)SA 

+(88.27D3A|-0.343)]-14.29U2+17.7Tv 

•n  [D,  (216.  71+159.  65AQ-D  2  (139.  92-281. 1A9)-D,  (87.  79-179.  5A,)1W 

H  I         U 

XYY 

[D 1 ( 8928 . 96+6639 . 34A , )+Dg ( 2375 . 84-4773 . 08A , )+D 3 ( 1490 . 67 

iYYu 

-2996. 97A2)]^  ,  P, (1 . 474+0. 77A! ).     D3 (0. 659-0. 572A7 )x 

**  T  "F  —  T  a 

YY  YY 


+  tDi  5. 41 {A, (216. 71+79. 83A1 ) }-D2 { 0 . 686+A2 (139, 
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•140.548A2  )}-D3A2  (87  .  79-88  .  27A2  )-0.  529-14  .  29U2  +17  .  7TX] 
IYY 

(3-24) 


From   (3-23)    and    (3-24) 


0    +   C26    -    C,Z   =   C36F   -    C^6A   +   U, 


-c6e  -  c7e  +  z  +  c5z  =  -c86F  -c9<sA  -  u2 


By  taking  the  Laplace  Transform  of  the  two  equations 
above  and  rearranging  into  a  matrix  form 


S(S+C2) 


-CS 


-(C6S+C7)      S(S+C5) 


-Ci 


-C. 


K 

+ 

'Uj  ' 

W 

-u2 

*                    * 

The  characteristic  equation 

A  =  S[S{S2+(C2+C5}-CiC6S-CiC7] 


wnere : 


C,  = 


[Dx(216.71  139. 65A1 )-D2 (139. 92-281. 1A, )-D3 (87 . 79-176 . 5A2 )] 

iYYu 


[D2 (8928. 96+6639. 84 Aj )+D2 (2375. 84-4773 . 08 A2 ) 


iYYu 


+D3 (1490. 67-2996. 97A2 )] 


P   -  [Di(l. 474+0. 77Ai)1 
XYY 


„      _    [D3(0. 659-0. 572A2)1 

YY 


Cs  = 


_  [5.162D!  +8.24D2+5.17D3 ] 


mU 


n  [(214.69Di-139.92D2-87.79D3 )  +U] 

C6  ~y 
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m 


«      _    [Di(0.0015Ai+0.036)] 

V--8   - — 


m 


After  inserting  the  value  of  parameters,  a  speed  of  36 
knots  was  chosen  to  be  used  in  this  simulation  to  insure 
that  cavitation  did  not  exist,  and  it  was  assumed  that  there 
was  no  external  disturbance.   The  transfer  function  matrix 
becomes : 


S(S+10.791) 


-0.063S 


-(76.33S+32.2)     S(S+7.61) 


Z 


0.083 


-1.04 


>  / 


-0.0394 


-1.215 


and  the  characteristic  equation  was 


A  =  S[S3+18.4S2+77.31S-2.03] 


and 


=  S(S-0.0261)(S+6.55)(S+11.876) 


0.083 

-0.063S 

6A=0 

-1.04 

S(S+7.61) 

0.083(S+6.82) 


(S+0.0261)(S+6.55)(S+11.876) 


6„ 


6A=0 


S(S+10.731)    0.083 


(76.33S+32.2)    1.04 


1.04(S-0.495)(S+5.19) 

S(S-0.0261)(S+6.55)(S+11.876) 
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e 
1 


6F=0 


•0.0394      -0.063S 


-1.215     S(S+7.61) 


-0.0394(S+9.55) 


(S-0.0261)(S+6.55)(S+11.876) 


6F=0 


S(S+10.791)    0.0394 


-(76.33S+32.2)   1.215 


=   1.215(S+0.079)(S+13.186) 

S(S-0.0261)(S+6.55)(S+11.876) 


W^  =  -Usin8+Vcos8sincj)+Wcos8cos4> 


V,  the  velocity  in  Y-direction  is  equal  to  zero,  and 
angle  0  is  assumed  small. 


Sin  6=6, 


Cos  6  =  1, 


w„  =  -u«e  +  w 


w 

"hcg  -  "ZE  "  "(S 


1-e, 


The  height  at  the  sensor, 


-h0  =  -h   +L„TT« 
S      eg   XH 


W  .  W  n     .  T 
~hS  =  -  S  +  S*6  +  LXH" 


-hS 


6r. 


6A=0 


+ 


u  _e_ 

s  6r 


6Z=0 


+  L 


e 


XH   6. 


6A=0 


6A=0 


-hS 


-_Z 
6. 


6F=0 


,  U   6 
+  S   67 


6F=0 


6F=0 


LXH   6A 


6F=0 
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The  four  transfer  functions  of  the  craft  are 


0.083(S+6.82) 


6A=0 


(S-0.0261)(S+6.55)(S+11.876) 


-hc 


6A=0 


6.68(S+0.72)(S+6.52) 
S(S-0.0261)(S+6.55)(S+11.876) 


6F=0 


-0.0394(S+9.55) 

(S-0.0261)(S+6.55)(S+11.876) 


-h. 


6F=0 


-1.464(S+2.72)(S+5.348) 
S(S-0.0261)(S+6.55)(S+11.876) 


The  characteristic  equation  shows  that  this  system  has 
one  pole  .in  the  right  half  plane,  and  from  the  Root  Locus 
plots  three  of  the  plant  transfer  functions  ha.ve  roots  in 
right  half  plane,  only  the  effect  of  forward  foil  to  the 
height  sensor  has  complex  roots  in  the  left  halp  plane  but 
very  close  to  the  imaginary  axis,  which  was  only  marginally 
stable.   From  this  point,  the  stability  of  the  craft  without 
automatic  control  system  was  questioned,  and  led  to  the  con- 
clusion that  the  craft  must  have  an  automatic  control  system. 

B.   ROOT  LOCUS  DESIGN  OF  FORWARD  FOIL  AUTOMATIC  CONTROL 

Considering  only  the  forward  foil  as  a  control  surface, 
the  aft  foil  was  fixed  at  the  angle  of  0.074269  radian  [Ref. 
1].   By  closing  the  height  loop  and  the  pitch  angle  loop  as 
in  Figure  14  and  increasing  gain  K  of  the  system  to  5,  from 
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Figure  10  the  roots  of  the  height  loop  moved  further  into 
the  left  half  plane  and  became  repeated  real  roots  at  -1.6. 
This  value  of  gain  K  pushes  the  root  of  the  pitch  angle  loop 
(Figure  11)  just  across  the  imaginary  axis  to  the  left  half 
plane.   The  system  becomes  stable.   The  response  to  the 
change  in  height  command  of  2.0  feet  was  shown  in  Figure  15. 
The  craft  went  to  steady  state  in  5.5  sec. 


?W  + 


'  — 


JL 


S 


e 


FIGURE  14 
BLOCK  DIAGRAM  SHOWING  FORWARD  FOIL  AUTOMATIC  CONTROL 

C.   ROOT  LOCUS  DESIGN  OF  AFT  FOIL  AUTOMATIC  CONTROL 

In  like  manner,  if  the  aft  foil  is  used  as  a  control 
surface  and  the  forward  foil  is  fixed  at  the  angle  of 
0.03711  radian,  the  system  is  unstable  by  itself,  because 
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FIGURE  15 
STEP  RESPONSE  OF  THE  CRAFT 
(FORWARD  FOIL  ONLY) 
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the  open  loop  roots  of  both  height  and  pitch  angle  were  in 
the  right  half  plane.   From  the  Root  Locus  plots  Figures  12, 
13,  positive  gain  will  have  the  roots  further  into  the  right 
half  plane.   After  closing  the  loops  and  using  gain  K  equal 
to  -47,  the  height  loop  had  complex  roots  at  -3±j3.   The 
right  half  plane  root  of  the  pitch  angle  loop  moved  toward 
the  left  half  plane,  across  the  imaginary  axis  to  the  point 
-0.2.   The  system  becomes  stable.   (See  Figure  16.) 
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—V 


FIGURE  16 
BLOCK  DIAGRAM  SHOWING  AFT  FOIL  AUTOMATIC  CONTROL 

Figure  17  shows  the  response  of  the  craft  to  the  change 
in  height  command  of  two  feet.   Because  the  flap  angle  was 
limited  at  ±5°,  and  the  gain  of  the  aft  foil  was  too  high, 
in  the  first  five  seconds  the  flap  exceeded  its  limit  and 
the  flap  angle  was  set  to  5°  which  had  the  same  result  as 
reducing  the  gain  of  the  system.   Then  in  the  first  five 


56 


,00 


12.  as 


id  oa 


23.  j: 


30.  o; 


X-Scale   6.0 
Y-Scale   1.0 


FIGURE  17 
STEP  RESPONSE  OF  THE  CRAFT 
(AFT.  FOIL  ONLY,  INPUT  =2.0) 
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seconds  the  craft  behaved  like  a  very  lightly  damped  system, 
but  after  five  seconds  the  .craft /became  more  stable  and 
reached  steady  state  in  eight  seconds.   Figure  18  shows  the 
response  of  the  craft  when  the  change  in  height  command  was 
0.5  foot.   In  the  first  1.0  sec.  the  flap  angle  still  ex- 
ceeded the  limit,  but  the  craft  went  to  steady  state  in  3 . 0 
seconds. 

Parts  B  and  C  have  shown  that  the  stability  of  the  craft 
in  pitch-heave  mode  can  be  achieved  by  using  either  forward 
or  aft  foil  as  a  control  surface  with  the  other  foil  fixed. 
But  to  make  the  craft  completely  stable  in  six  degrees  of 
freedom  both  forward  and  aft  foil  have  to  be  used  as  control 
surfaces.   The  dynamics  of  the  system  will  produce  a  cross 
coupling  .effect  between  the  two  inputs,  the  systems  now  have 
multiple  inputs  and  multiple  outputs,  which  is  a  classical 
multivariable  system  and  will  be  studied  in  the  next  chapter 
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FIGURE  18 
STEP  RESPONSE  OF  THE  CRAFT 
(AFT.  FOIL  ONLY,  INPUT  -  0.5) 
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V.   MULTIVARIABLE  SYSTEM 

A.   BASIC  CONCEPT 

A  system  which  has  multiple  inputs-multiple  outputs  is 
much  more  complicated  to  design  than  a  single  input-single 
output  system,  because  a  change  in  one  input  may  have  an  ef- 
fect not  only  on  the  corresponding  output  but  also  on  the 
other  outputs.   In  designing  this  kind  of  system,  special 
techniques  are  needed.   Ref .  9  has  all  of  the  theories  and 
procedures  to  work  on  this  kind  of  problem,  and  this  paper 
will  follow  the  procedures  in  that  reference. 

The  idea  was  that,  the  plant  has  a  nxm  transfer  function 
matrix  and  was  compensated  by  a  iPn  diagonal  matrix  compen- 
sator with  unity  feedback.   The  configuration  of  the  compen- 
sated system  is  shown  in  Figure  19. 

Consider  the  hydrofoil  craft  in  pitch-heave  made  only, 
and  not  coupled  to  the  other  modes.  The  craft  will  have  a 
2x2  transfer  function  matrix. 


G   = 
P 


g 


Pn 


'2  1 


Pi  2 


P 


2  2 


(4-1) 


If  the  diagonal  compensator  matrix 


G   = 
c 


e. 


■i  i 


o 


0 


C2  2 


(4-2) 
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is  used,  the  system  configuration  becomes 
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FIGURE   19 

BLOCK  DIAGRAM  SHOWING  A  MULTIVARIABLE  SYSTEM 
WITH  DIAGONAL  MATRIX  CASCADE  COMPENSATOR 


G  G   = 
P  c 


gPnSCn 


P2  1   Cj.  1 


gPl2gC22 


P2  2   C2  2 


(4-3) 


1.   The  Characteristic  Equation 

The  stability  of  a  multivariable  system  is  determined 
by  the  zeros  of  Ni(s)  and  N(s),  which  are  defined  as  follows: 


f?1^?  A  det  [I+G  (s)G  (s)] 
Di  (s)  =         pv  J    cv  y 


N(s)  A 


Ac(s)-Ap(s) 

Di(s) 


(4-4) 
(4-5) 


and  A  (s),  A  (s)  represent  the  characteristic  polynomials  of 
c      p 

the  rational  transfer  function  matrices  G  (s)  and  G  (s). 

cv        p   ' ' 

respectively. 
If 

1)   Cancellations  are  selected  systematically,  as  in 
Ref.  9, 
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2)   Poles  of  G   are  carefully  selected, 
c 

then  (4-4)  alone  determines  the  stability  of  the  system. 
Thus 

det(I+G  G  )  =  0 
p  c 

will  be  referred  to  as  the  characteristic  equation  of  multi- 
variable  system. 

2.   Connection  Between  Single-Loop  and  Multivariable 
Cases  for  2x2  Plant  G  ,  and  2x2  Diagonal  G 


det(I+G  G    )    =   1+g        g        +g        g        +(det-G    )g        g 

P     C  Pi  1       cl  1  P2  2      c2  2  P         cl  1       c2  2 

g        +(det-G    )g 

=       d+g„  g„  )       1+     -J^-ib rr  M        g 


Pn  Cii        l+gn  Sr  &c22 

Pi  i  ci  i 


where 


(4-6) 


=  (1+g    gc   )d+G   g    )  (4-7) 

Pn  cn      eq  c2  2 


det-G 
g   +(det.G  )g      gp22  1+   gn    gcn 

eq     1+g   g  1+g   g 

Pi  1  ci 1  Pi  1  ci 1 


The  roots  of  the  characteristic  equation  are  the 

zeros  of  the  rational  function  (1+g   g    )(1+G   g    ). 

Pi  1  Ci 1      eq  c2  2 

From  (4-6),  can  be  seen  that  the  roots  of 

det«G 

1  +  E  g    =0  (4-8) 

g„      Cn  v 

P2  2 


and 


1  +  gn  Sn        =0  (4-9) 

Pi  1  Ci  1 
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constitute  part  of  the  zeros  and  poles  of  G   ,  respectively 

det'G 
and  2  and  g    are  known  functions. 

gP22       Pl1 

Therefore,  any  single-loop  method  can  be  used  in  de- 
signing g    to  place  the  roots  of  (4-8)  and  (4-9)  at  desir- 

Ci  l 

able  locations,  and  also  has  toform  a  reasonably  good 

pole-zero  pattern  for  G   .   After  g    is  designed,  all 

*  *  eq  cii 

poles  and  zeros  of  G   are  known  and  g    can  be  designed 

*  eq  &c22 

from 


1+G   g    =0 
eq&c22 


(4-10) 


B.   ROOT  LOCUS  DESIGN  FOR  PITCH-HEAVE  MODE 
From  the  craft  transfer  function  matrix 


G  = 
P 


-h. 


I  6F 


-h. 


'A 


g 


Pa  i 


6.68(S+0.72)(5+6.52) 

S(S-0.0261)(S+6.55)(S+11.876) 


g 


Pl2 


A 


-1.464(S+2.72)(S+5.348) 
S(S-0.0261)(S+6.55)(S+11.876) 


g 


P2  1 


0.083(S+6.82) 


(S-0.0261)(S+6.55)(S+11.876) 


g 


P22 


-0.0394(S+9.55) 


(S-0.0261)(S+6.55)(S+11.876) 
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det-G^  =  g^  g^  -g   g^ 

P    Pi  1  P22   P12  P21 

_   -0.141(S-0.02325)(S+6.51)(S+12.02) 
S(S-0.0261)2(S+6.55)2(S+11.876)2 

det ' Gp    3.579(S-0.02325)(S+6.51)(S+12.02) 

g      ~   S(S-0.0261)(S+6.55)(S+11.876) 
P22 

The  design  procedure  starting  from  g    ,  by  plotting  the 

Ci  1 

root  locus  of  equations  (4-8)  and  (4-9)  as  in  Figures  20  and 
21,  from  Figure  20  can  be  seen  that,  as  the  gain  was  in- 
creased the  root  in  the  right  half  plane  moved  toward  the 
left  half  plane,  met  pole  at  the  origin,  broke  away  and 
circled  around  the  zero  at  -0.72,  with  the  value  of  the  gain 
equal  to  5.   Two  dominant  roots  were  at  -1.25  and  -2.2,  and 
other  two  roots  at  -6.75  and  -11.45  and  the  effect  of  this 
gain  pushes  the  right  half  plane  root  of  equation  (4-8)  to 
the  left  and  was  cancelled  by  the  zero  at  +0.02325  as  in 
Figure  21,  the  pole  at  the  origin  moved  further  into  the 
left  half  plane  to  the  point  -2.3  and  the  other  roots  were 
in  the  left  half  plane  at  -6.8  ±  j0.8  and  -11.99.   Therefore, 

by  using  the  gain  equal  to  5  as  g    ,  the  height  loop  was 

Ci  1 

stabilized  and,  this  also  gave  a  good  pole-zero  pattern  for 

G   . 
eq 

Equations  (6-30)  -  (6-35)  of  Ref .  9  show  the  procedures 

to  find  poles  and  zeros  of  G   ,  which  were  listed  below 

eq 

Gi  =  l+gn   gp  (4-11) 

Pi  1  Ci  1 

G2  =  gn   +(det-G  )g  (4-12) 

P2  2         P   Ci  1 
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Geq      =§f  <4"13> 

V  ■zG1+por<pQ»npol)  (4"14) 

6  (J 

V   =ZG2+PGr(PG,   PG,'  <4"15> 

P„    =  P            =  P     +P  *    (4-16) 

Gl      gpngcn        gpn   gCll 


Zr    =  Zi+g    g  (4-17) 

<->i  Pi  1  Cn 


det  •  G 

Zr    =  Zi  +  £  g  (4-18) 

G2  g„     en 

P2  2 

And  from  the  system  transfer  function  matrix 

PQ      =   P-  P  =   P  +P  =    {0,0. 0261,-6. 55, -11. 876}+P 

Piicii  Piicii  cii 


ZG      =    Zl+gt)      gc 
ui  Pi l    ci l 

P„      =  P  +P(det«G    )g        -[P       fl      P(det«G    )g        ] 

G2  gn  v  P;5cu         g  v  p;5C!i 

p2  2  P2  2 

=    {0.0261,-6. 55,-11. 876} +{0, (0.0261)2 , (-6 . 55 )2 , (-11 . 876  )2 } 

+P  -[{0.0261,-6.55,-11.876}D   { 0 , ( 0. 0261 )2 , (-6. 55)2 , 

gcn 

(-11.876)2}+gr>       ] 
ci  i 

=    {0, (0.0261)3 ,(-6.55)3 ,(-11.876)3    +P  -[{(0.0261), 

gCn 

(-6.55), (-11.876)}] 
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=  {0, (0.0261)2 , (-6.55)2 , (-11.876)2}+P 

gcn 

det*G 
Zn      =  Zi+  £  g 

P2  2 

P^  fl  P^   =  {0,0. 0261,-6. 55, -11. 876}+P      =  P„ 
Oi    b2  g        G 

Ci  i       l 


Pr    =  Zi+g„   g„   +  {0,(0.0261)2,(-6.55)2,(-11.876)2}+Po. 
Geq       Pl1  Cl1  gCll 

-{0,0. 0261,-6. 55,-11. 376}-P' 

gcn 

=  Zi+g^   g„   +{(0.0261), (-6.55), (-11.876)} 

Pl 1  Ci i 

det-G 

Z     =  Z     =  Zx+  — ^  g^ 

eq     z         &p22     X1 


G   has  three  known  poles  at  0.0261,  -6.55,  -11.876;  the 
eq 

other  poles  were  the  roots  of  equation  (4-9) 

l+gn   g_    =0 
Pi i  ci i 


The  zeros  of  G   were  the  roots  of  equation  (4-8) 
eq  ^ 

det-G 

1+  E  g    =0 

g„     en 

P2  2 


and  therefore 


K   (S+2.3)(S+6.8±j0.8)(S+11.99) 


eq    (S-0.0261)(S+1.25)(S+2.2)(S+6.55)(S+6.75)(S+11.45)(S+11.876) 
and  from 
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g^   +(det-G  )g^ 

G     _    P2  2  P    Cl  1 

eq  "   l+g„   g„ 

Pi  1  Ci  l 

0.0394(S  +  9.55) 0. 144  (S-0.  02325)  (S+6  .  51)  (S+12  .  02) 

(S-0.0261)(S+6.55)(S+11.876)~  S(S-0.0261)2 (S+6.55)2 (S+11.876)2 

1     6.68(S+0.72)(S+6.52) 

S( S-0. 0261) (S+6. 55) (S+11.876) 

-0.0394(S+Z1 )(S+Z2 ) 

(S+Pi)(S+P2) 


K   =  -0.0394 
eq  (4-19) 

r    =  -0.0394(S+2.3)(S+6.8±j0.8)(S+11.99) 

eq    (S-0.0261)(S+1.25)(S+2.2)(S+6.55)(S+6.75)(S+11.65)(S+11.876) 

Equations  (4-10)  and  (4-19)  were  the  equations  used  in 
designing  the  compensator  g 

P2  2 

Figure  22  is  the  root  locus  plot  for  designing  g 

P2  2 

The  positive  gain  will  push  the  pole  at  0.0261  further  into 
the  right  half  plane;  therefore  positive  gain  cannot  be  used 
to  compensate  this  system.   Negative  gain  which  can  make  the 
gain  of  the  system  equal  to  3.881  will  set  dominant  roots  at 
-0.7±j0.2,  which  is  a  pretty  good  position.   Therefore, 
using  the  gain  of  the  system  equal  to  3.881,  we  can  make 
system  stable  and  the  pitch  angle  loop  has  good  damping. 
The  compensator  g    thus  becomes  a  pure  gain  of  value 

P2  2 

K2  K    =  3.881 
eq 

K2  =  -o!o396  =  ~98 
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The  configuration  of  the  compensated  system  is  shown  in 
Figure  23. 
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FIGURE  23 
BLOCK  DIAGRAM  SHOWING  THE  COMPENSATED  SYSTEM 

Referring  to  Figure  20,  the  dominant  roots  were  at  -1.25 
and  -2.2  which  gave  co   =  1.658,  and  ;;  =  1.04.   Figure  2-1 
shows  the  response  of  the  craft  when  the  height  above  the 
water  has  changed  from  9.5  feet  to  11.5  feet,  and  Figure  25, 
when  it  changes  from  11.5  feet  to  9.5  feet.   The  response  of 
the  craft  was  the  same  in  both  cases  and  followed  very 
closely  to  the  response  of  the  best  second  order  approxima- 
tion. 

Figure  26  is  the  response  of  the  pitch  angle  from  2°  bow 
up  back  to  0°  reference  and  Figure  27  is  the  response  from 
2°  stern  up  back  to  0°  reference.   The  response  in  both 
cases  was  the  same  but  faster  than  the  response  of  the  best 
second  order  approximation. 
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FIGURE  24 
STEP  RESPONSE  OF  HEIGHT  (2  FT.  UP) 
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FIGURE  25 
STEP  RESPONSE  OF  HEIGHT  (2  FT.  DOWN) 
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FIGURE  26 
STEP  RESPONSE  OF  PITCH  ANGLE  (BOW  UP) 
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FIGURE  27 
STEP  RESPONSE  OF  PITCH  ANGLE  (STERN  UP) 
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C.   ROOT  LOCUS  DESIGN  FOR  ROLL  MODE 

After  neglecting  the  small  terms,  equation  (2-4)  becomes 

XX 

And  from  equation  (2-41) 

L  "  FZPFLYPF+FZSFLYSF~FYPSLZPS~FYSSLZSS_FYCSLZCS 

Assume  that  the  roll  mode  was  decoupled  from  the  other 
modes 

L   =  FZPFLYPF+FZSFLYSF 

=  [-D2(0.0028A2+0.018)6p]LYpF+[-D2(0.0028A2+0.018)6s]LySF 

6P  =  "6S 

L  =  -29959. 46s+(-29959.4<5g) 
=  -59918. 86g 

P   =  cj>   =  -O.lfig 

Taking  the  Laplace  Transform,  the  transfer  function  is 
simply 

%•-*$  (4-20) 

Figure  28  is  the  root  locus  plot  after  compensation 
using  a  lead  network  which  has  zero  at  -1.5  and  pole  at  -15. 
With  the  gain  of  -82,  the  roots  were  at  -3.5±j0.75,  which 
gave  the  value  of  £  =  0.96  and  to   =  3.6.   Figure  29  is  the 
step  response  roll  from  2°  port  up  back  to  0°  reference  and 
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Figure  30  from  2°  starboard  up  to  0°  reference.  The  same 
response  is  obtained  in  both  cases.  The  block  diagram  of 
the  automatic  control  of  this  system  is  shown  in  Figure  31 
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FIGURE  29 
STEP  RESPONSE  OF  ROLL  ANGLE  (PORT  UP) 
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FIGURE  30 
STEP  RESPONSE  OF  ROLL  ANGLE  (STARBOARD  UP) 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  transfer  functions  and  root  locus  study  has  shown 
that  this  system  is  unstable  by  itself  because  of  the  open 
loop  pole  in  the  right  half  plane.   From  studying  the  effect 
of  each  foil  on  the  system,  it  has  been  shown  that  stabiliza- 
tion can  be  achieved  by  simply  closing  each  loop  and  using 
the  appropriate  gain.   However,  that  will  stabilize  only  the 
effect  of  one  foil  when  the  other  is  fixed. 

To  stabilize  the  systems  when  both  forward  and  aft  foils 
are  used,  the  total  system  has  been  treated  as  a  multivari- 
able  system.   To  reduce  the  work  in  designing  the  automatic 
control  system,  some  assumptions  had  to  be  made,  such  as  no 
coupling  between  the  pitch-heave  modes,  roll  mode  and  yaw 
mode.   Thus,  the  system  will  have  only  a  2^2  transfer  func- 
tion matrix  for  the  pitch-heave  mode,  and  single  loop  for 
the  roll  mode.   The  step  response  study  has  shown  that  the 
responses  of  the  compensated  system  followed  closely  the  re- 
sponse of  the  best  second  order  approximation.   The  system 
was  made  stable  in  both  pitch-heave  and  roll  mode.   However 
this  result  cannot  insure  that  the  system  is  completely 
stable  in  six  degrees  of  freedom. 

Recommendations  for  further  study  are: 

1)  Study  of  yaw  motion  and  the  craft  turning  character- 
istic 

2)  Wind  effects 
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3)  Produce  a  simulation  for  the  transition  period  be- 
tween hullborne  and  foilborne 

4)  Take  the  dynamic  characteristic  of  the  sensors  into 
account 

5)  Study  the  response  of  the  craft  in  a  sinusoidal  sea 
and  in  a  random  sea. 
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APPENDIX  A 


SUMMARY   OF   COMPLETE    EQUATIONS 


EQUATIONS    FOR   ACCELERATION    IN   BODY   AXIS 

U   =  -F   +RV-QW  (A-l) 

m  X 

V  =  -F..+PW-RU  (A- 2) 

ra  Y 

W  =  ^F   +QU-PV  (A- 3) 

P   =  T^-[L-QR(IZZ-IYY)+(R+QP)IXZ]  (A-4) 

XX 

R  =  I^-[N-PQ(IYY-IXX)-(QR-P)IXZ]  (A-5) 


VELOCITY    IN   BODY   AXES 

U   =   /   Udt  (A-6) 

V   =   /   Vdt  (A~7) 

W   =    /    Wdt  (A-8) 

P   =   /   Pdt  (A-9) 

Q  =    /   Qdt  (A-10) 

R   =    /   Rdt  (A-ll) 

TRANSFORMATION   OF   VELOCITIES    FROM   BODY   TO  EARTH  AXES 

Uv=Ucos0cos*jj+V(cosiJ;sin6sin(j)-siniJ;cos(J)  )+W(cosiJjsin0cos(}>+siniJjsin(j) ) 

*  (A-12) 

Vp=Ucos0sim(;+V(cosiJjcos(f)+sirn|jsin6sin<{))+W(siniJjsin0cos(t)-cosiJjsin4)) 

B  (A-13) 

W  =-Usin0+Vcos0sinc|)+Wcos0cos(J)  (A-14) 
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<J>    =   P+^sin9  (A-15) 

0    =  Qcos<f>-Rsin<{>  (A-16) 

i>   =    (Qsin4>+Rcos<t>)cose  +  (<i>-P)sin9  (A-17) 

POSITION  IN  EARTH  AXES 

X£   =   /   UEdt  (A-18) 

YE   =   /   VEdt  (A-19) 

Z£   =    /   Y;£dt  (A-20) 

<j>      =   /   4>dt  (A-21) 

Q      =   f    8dt  (A-22) 

ty      =    /    Tpdt  (A-23) 

EXPANSION   OF   FORCE   AND   MOMENT   EQUATIONS 

FX   =    E    FXIF    +   E    FXIS    +   mgX+TX  (A"24) 

Fy   =   £    Fyig   +   mgY  (A-25) 

Fz   =   Z    FZIF   +  mgz  (A-26) 

L      =    (FzLY)pF+(FzLY)SF-(FYLz)ps-(FYLz)ss-(FyLz)cs    (A-27) 

M   =  Z    (FZLX)IF+E  (FXLZ)IF+E  (FXLZ)IS+TXLZT  (A-28) 

N   =  Z    (FyLx)Ig  (A-29) 

EXPANSION  OF  GRAVITY  TERMS 

mg   =  -mg  sine  (A-30) 

mgy  =  mg  cosesintj)  (A-31) 

mg„  =  mg  cosesint})  (A-32) 
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FOIL  VELOCITY  COMPONENTS  IN  BODY  AXES 


Center  Foil 


Uc  =  U+LZCF-Q  (A- 33) 

VC  =  V"LZCF- P+LXCF* Q  <A-34> 

Wc  =  W-Lxcp-Q  (A-35) 


Port  Foil 


UP  =  U+Lzpf'Q-LYPF*R  (A"36) 

VP  =  V-LZPF'P+LXPF,R  (A_37) 

Wp  =  W+LYPF*P-LXPF,Q  (A"38) 


Starboard  Foil: 

US  ="  U+LZSF*Q-LYSF*R  (A~39) 

VS  =  V-LZSf'P+LXSF'R  (A"4°) 

WS  =  W+LYSF*P-LXSF'Q  <A"41) 


Mid  Foil: 

.  UM  =  U+LXMF*Q  <A"42) 

VM  =  V-LZMF'P+LXMF'R  (A~43) 

WM  =  W"LXMF-Q  (A"44> 
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TRANSFORMATION   OF   WATER  PARTICLE   ORBITAL   VELOCITY   FROM   EARTH 
TO   BODY   AXES 

UWI    =   uEWicos^cos9+vEWisin^cos0-WEWISin9  (A-45) 

VWT   =  UTrWT(cosi|;sin0sin(f)-sin4)COS({)  )+VWEI(siniJjsinOsin(f)+cosiJjcos(})) 

+WEWIcos9sincj)  (A-46) 

WWT   =  U_WT(cos^sin9cos<j)+sin^sin<J)  )+VWEI(sin^sin6cos(J)-cos^sin4)) 

+WEWIcos6cos(()  (A-47) 

Each  of  these  equations  must  be  repeated  for  each  foil. 

RELATIVE  VELOCITY  COMPONENTS 

URI  -  Uj  -  UWJ  (A-48) 

VB1"-  Vj  -  VWI  (A-49) 

W    =  W   -  W  (A-50) 

RI     I     WI  V 

One  set  of  equations  for  each  foil. 

ANGLES  OF  ATTACK  AND  SIDE  SLIP 

Oj  =  arctan  WRI/DRI  +  oI_FIXED  (A-51)* 

VRT 
6T  =  arcsin      K1   (A-52)* 

1  /U2  +V2  +W2 

Y    RI   RI   RI 

*  One  set  of  equations  for  each  foil.   Angles  of  attack  and 
side  slip  are  calculated  in  radians. 
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TOTAL  RELATIVE  VELOCITY  AT  A  PARTICULAR  FOIL  OR  STRUT 


V  =  /U^+.V^+W^  (A-53) 


HYDRODYNAMIC  FORCES  IN  WATER  AXES 

(Each  foil  equation  adn  strut  equation  must  be  repeated  for 

foil  and  strut,  respectively) 

Foil  Lift: 

FL   =  -JpV2AFCF  (A- 54) 

Foil  Drag: 

FDF  =  *PV2AFCF  (A_55) 

Strut  Drag: 

FDS-  *PV2ASCDS  (A-56) 

i 
Strut  Side  Force: 

Fg   =  ipV2AgCs  (A-57) 

TRANSFORMATION  OF  FOIL  FORCES  FROM  WATER  AXES  TO  BODY  AXES 

FXIF  =  -FIDFcosacosP+FL];Fsina  (A-58) 

FYIF  =  -FDIFSinB  (A_59) 

FZIF  =  -FDIFsinacos3-FLIFcosa  (A-60) 
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APPENDIX  B 


PLOTS  OF  HYDRODYNAMIC  COEFFICIENTS 


Fwd.  Foil  Lift  Coefficient  vs.  Angle  of  Attack 


-8-do  -wo  4-00 

X-Scale  4,0  degrees/inch 

Y-Scale  0.4  Unit/inch 


8  OO 


iZ-oet 
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Mid  Foil  Lift  Coefficient  vs.  Angle  of  Attack 


400 


BOO 


X-Scalo     4,0     Degrees/inch 
Y-Soal©     0o4     Unit/inch 
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Aft.  Foil  Lift  Coefficient  vs.  Angle  of  Attack 


r 

-800?  -A-OO 

X-Scale  4.o  Degrees/inch 


A-OO 


B-OO 


lzoo' 


Y-Scale  0.4  Unit/inch 
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Forward  Foil  Drag  Coefficient  vs.  Angle  of  Attack 


o 

o 

Ha 


-* 


N 


o 

u-o 


*^0 


1  Flap=-5.0 

2  Flap=  0.0 

3  Plap=+5.0 


1 1 1 — 

8-00  -4-0O  4--00 

X-Scale     4.0      Degrees/inch 
Y~Scale     0.02  Unit/inch 


B-OO 


izoo 
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//  EXEC  DSL,PARM.C=,NOSOURCE,NOMAP' 

//DSL. INPUT  DD  * 

D  CCMMON    MATRIX(539) 

CONTRL    FINTIM=10.0,DELT=0.025,DELS=0.05 

INTEG    RKSFX 

INTGER     IFLAG, I F LAG1 , NPLOT ,NUM, MCUR 

CONST    NPLOT=3 ,NCUR=1,IFLAG=0 

CONST    WEIGHT=262527,     IXX=601200,     IYY=4286000, . . . 

IZZ=41 7400G, G  =  32. 2 

CONST    ACFP=32.7,ACFS=32.7,APF=36.8,ASF=36.3,AHF=76.0,    ... 

PI=3. 141593 

CONST    LXCF  =  41 .59 , LXCS  =  41 . 59, LXPF  =  -1 6 . 9  8 , LXPS  =  - 16 . 93 ,  .  .  . 

LXSF=-lo.98,LXSS=-16.9  8,LXMF=-16.98,LYCP=-5.0, LYCS=5.0,. . 

LYPF=-12.42,LYPS=-8.0,LYSF=12.42,LYSS=8.0,LZCF=15.5,... 

LZCS=13.25,LZPF=17.14,LZPS=14.0,LZSF=17.14,LZSS=14.0,... 

LZMF=17.14,LXT=-16.9,LZT=17.7,LXH=6  8.0 

CONST    Kl=5.0,K2=-98.0,K3=-82.0 

STORAG     IC(2I tNU(l) f0EN(3) 

TABLE    IC(l-2)=2*0.0,NU(l)=3.99iDEN{l-3)=1.0,4.0,3.99 

PARAM    UZERO=36.0 

PREPAR     .1,FX,FY,FZ  ,L, M , N , U, V , W, Z E ,P , Q , R , THETA     ,  UDOT  ,V^DOT  ,  . 

QOOTt ALFCPi ALFCSf ZElf PHi , PH 1 1 , THE TA3 , X , Y 

GRAPH    TIME,THETA3,X,Y 

PRPLOT 

LABEL  STEP  RESPONSE  OF  PITCH  ANGLE 

PRINT  0.5, TIME, Z E, PHI ,P , PDOT , THETA, EL VP , ELVS , F LAP , ARS , ARP 

* 

INITIAL 
* 

DENSE=2.0 

UINIT=U ZERO* 2 000. 0*3. 0*1. 0/3600.0 

U=UINIT 

TX=3960.0 

A LM=0. 074269 

ALS=0. 074269 

ALP=0. 074269 

ALF=0. 037111 

PREF=0.0 

TREF=0.0 

THETA2=-2.0*PI/180.0 

PHI2=0.0 

MASS=WEIGHT/G 

CONV=PI/180.0 

FLAP=0.0 

ELV=0.0 

ARP=0.0 

ARS=0.0 

HREF=11  .5 

HREF2=-11.5 

FLAG=0.0 

IFLAG1=0 

EPS=25.0 

DERIVATIVE 

NOSORT 

*      ANGLES  OF  ATTACK  AND  SIDE  SLIP 


* 


VRC=V-LZCF*P+LXCF*R 

VRS=V-|_ZSF*P+LXSF*R 

VRM=V 

VRP=V-LZPF*P+LXPF*R 

WRP=W+LYPF*P-LXPF*Q 

WRS=W+LYSF*P-LXSF*Q 

WRM=W-LXMF*Q 

WRC=W-LXCF*Q 

ALFC=WRC/U+ALF 

ALFAP=WRP/U+ALP 

ALFAS=WRS/U+ALS 

ALFAM=WRM/Uf ALM 

6ETAC=VRC/U 
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BETAP=VRP/U 

BETAS=VRS/U 
* 

*      CALCULATION  OF  EULER  ANGLES 
* 


* 
SORT 


PHID1=P 

THETD1=Q 

PSID1=R 

THETA=INTGRL(THETA2,THETD1 J 

PHI=INTGRL(PHI2, PHID1) 

PSI=INTGRL(0.0,PSID1) 

WE=-U*SIN(THETA) +V*COS ( THET A )*S IN ( PHI ) +W*COS< THET A ) 

*COS(  PHI) 

ZE=INTGRL(HREF2,WE) 

ZE1=ZE-HREF2 
* 
NOSORT 

*  SUBMERGENCE    TERMS 
* 

SECP=ZE-LXCF*THETA+LYCP*SIN(PHI ) +LZCF*C OS ( PHI ) 
SECS=ZE-LXCF*THETA  +  LYCS*SII\I(PHI  )  +LZCF*CGS  (  PHI  ) 
SEPF=ZE-LXPF*THETA+LYPF*SIN(PHI ) +LZPF*COS ( PHI ) 
SESF  =  ZE-LXSF*THETA+LYSF*SIN<PHI  )  +  LZi>F*COS  (  PHI  ) 
SEMF=ZE-LXMF*THETA+LZMF*COS(PHI  ) 
* 

*  HYDROOYNAMIC  COEFFICIENTS 
* 

ELVl=ELV*CONV 
ELVD=LI  MIT(-5.0,5.0,ELV1) 
FLAD=FLAP*CONV 
FDEG=LIMIT(-5.0,5.0,FLAD) 
ARPD=ARP*CONV 
ARP1=LIMIT(-5.0,5.0,ARPD) 
ARSD=ARS*CONV 
ARS1=LIMIT(-5.0,5.0,ARSD) 
CALL    INTERP(ALFC,ELVD,C0EFC,2) 
CALL    INTERP(ALFAP,ARP1 ,COEFPX,3  ) 
CALL    INTERP(ALFAS,ARS1,C0EFSX,3) 
CALL    INTERP(ALFAM,FLAP,COEFMX,l ) 
CALL    INTRPKBETACf  SECFtCFDSt2i  THETA,THETA1) 
CALL    IMTRPK  BET AP  ,  SEPF  ,CPDS , 1 , THETA , THE TA1 ) 
CALL    INTRPK  BET AS, SESF,CSOS,l, THETA, THET Al ) 
* 

*  FORCES    IN    THE    WATER    AXES 

FGRAX=-MASS*G*S IN (THETA) 

FGRAY=MASS*G*SIN(PHI) 

FGRAZ=MASS*G*COS( PHI ) 

CCS=-2. 15*BETAC 

ACCS=CCS*4.25*SECF 

FSC=0.5*DENSE*(U**2)*ACCS 

FDSC=0. 5*OENSE*( U**2) *CFOS 

PCS=-2. 15*BETAP 

AFCS=PCS*4.25*SEPF 

FPS=0.5*DENSE*(U**2)*APCS 

FDPS=0. 5*DENSE*(U**2) *CPDS 

SCS=-2.  15*BETAS 

ASCS=SCS*4.25*SESF 

FSS=0.5*CENSE*(U**2)*ASCS 

FDSS=0»  5*DEIMSE*(  U**2)  *CSDS 

CLCF=5.15#ALFC+0.036*ELV+0. 13 

FLC=0.5*DENSE*(U**2)*ACF*CLCF 

FDC=0.5*DENSE*(U**2)*ACF*C0EFC 

CLPF=4. 10*ALFAP+0.018*ARP 

FLP=0.5*DENSE*(U**2)*APF*CLPF 

FDP=0.5*DENSE*(U**2)*APF*COEFPX 

CLSF=4. 10*ALFAS+0.018*ARS 

FLS=0.5*DENSE*(U**2)*ASF*CLSF 

FDS=0.5*DENSE*(U**2)#ASF*COEFSX 
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CLMF=5. 15*ALFAM+0.036*FLAP 
FLM=0. 5*DENSE*(U**2)*AMF*CLMF 
FDM=0.5*DENSE*(U**2)*AMF*C0EFMX 
FLCZ=0. 5*DENSE*(U**2)*ACF*CLCF 
FDCZ=0.5*DENSE*(U**2)  *ACF*CQEFC 
FLPZ=0. 5*DENSE*(U**2)*APF*CLPF 
FDPZ=0. 5*DENSE*(U**2) *APF*COEFPX 
FLSZ=0.5*DENSE*(U**2)*ASF#CLSF 
FOSZ=0. 5*DENSE*(U**2)*ASF*CGEFSX 
FLMZ=0. 5*DENSE*<  U**2) *AMF*CLMF 
FDMZ  =  0. 5*DENSE*(U**2) *AMF*COEFMX 
* 

*  TRANSFORMATION    TO    BODY    AXES 

FXCF  =  -FDOALFOFLC 

FXPF=-FDP+ALFAP*FLP 

FXSF=-FDS*ALFAS*FLS 

FXMF=-FDM+ALFAM*FLM 

FXCS=-FDSC-BETAC*FSC 

FXPS=-FDPS-BETAP*FPS 

FXSS=-FDSS-BETAS*FSS 

FYCS=-BETAC*FDSC+FSC 

FYPS=-6ETAP*FDPS+FPS 

FYSS=-BETAS*FDSS*-FSS 

FZCF=^FDC*ALFC-FLCZ 

FZPF=-FDPZ*ALFAP-FLPZ 

FZSF=-FDSZ*ALFAS-FLSZ 

FZMF=-FDMZ*ALFAM-FLMZ 

LZCS=LZCF-0.4*SECF 

LZPS=LZPF-0.4*SEPF 

LZSS=LZSF-0.4*SESF 

FXONE=FXCF+FXPF*FXSF+FXMF+FXCS+FXPS+FXSS 

*  SUMMATION  OF  FORCES  AND  MOMENTS 
* 

FX=FXONE+FGRAX+TX 

UDOT=( 1 ./MASS)*FX-Q*W 
IF(  IFLAG1.EQ.0)  UDOF=0„0 
FY=FYCS+FYPS+FYSS+FGRAY 
FZ=FZCF+FZPF+FZSF+FZMF+FGRAZ 

L=FZPF*LYPF+FZSF*LYSF-FYPS*LZPS-FYSS*LZSS-FYCS*LZCS 
MFWD=FZCF*LXCF 

M=-MFWD-FZPF*LXPF-FZSF*LXSF-FZMF*LXMF  +  . .  . 
FXCF*LZCF+FXPF*LZPF+FXSF*LZSF+FXMF*LZMF+FXCS.  .. 
*LZCS+FXPS*LZPS*FXSS*LZSS+TX*LZT 
N=FYCS*LXCS+FYPS*LXPS+FYSS*LXSS 
* 

*  VELOCITY    INTEGRALS 

SORT 
* 

VDOT=(l ,/MASS)*FY-R*U 

WDOT=( 1 ./MASS)*FZ+Q*U 

PDOT=(1.07IXX)*L 

QDOT=(l .0/IYY)*M 

RDOT=(l  .0/IZZ)*N 

U=INTGRL(UINIT,UDOT) 

V=INTGRL(0.0,VDOT) 

W=INTGRL(0.0,WDOT) 

P=INTGRL(0.0,PDOT) 

Q=INTGRL(0.0,QDOT) 

R=INTGRL(0.0,RDOT) 
* 

*  AUTOMATIC    CONTROL    BLOCK 

jjs  if. 

HACT=-ZE+LXH#THETA 

DELH=HREF-HACT 

DETH=TREF-THETA 

DEPH=PREFxPHI 

ELV=K1*DELH 

FLAP=K2*DETH 
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ARS1=LEDLAG(0.0,0.66  7,0.066  7,DEPH) 

ARS=K3*ARS1 

ARP=-ARS 

THETA3=(THETA-THETA2)*180.0/PI 

PHI1=(PHI-PHI2)*180.0/PI 

X1=STEP(0.0) 

Y=TRNFR(0,2, IC,NU,DEN,X) 
* 

DYNAMIC 
* 

SAMPLE 

CALL    DRWG(1,1,TIMF.,THETA3) 
CALL    DRWGd  ,2, TIME, X) 
CALL    DRWGi 1,3, TIME, YJ 

TERMINAL 

CALL  ENDRW(NPLOT) 
END 
STOP 
FORTRAN 
C 
C 

C  SUBROUTINE  SETUP 

C      THIS  SUBROUTINE  FORMS  THE  DATA  POINTS  COMPRISING 
C      THE  CURVES  OF  FOIL  AND  STRUT  DRAG  COEFFICIENTS 
C      AND  SIDE  SLIP  COEFFICIENTS  INTO  MATP.ACES. 
C 

c 

SUBROUTINE  SETUP(FLAG) 

COMMON  COEFMX,COEFCX,COEFOX,CFDSl,CADSl 
DIMENSION  COEFCX(3,33) , COEFOX ( 3 , 33 ) , COEFMX ( 3 , 33 ) ; 
lCFDSnilill)  fCADSlIll  tll)f  INFO(  8) 
REAL  INFO 
IC  =  0 
IFLAG=0 
IFLAG1=0 
tfATRIX=l 
J=0 
K=l 
ICOUNT=0 

2  READ(3,100,END=50)  INFO 
100  FORMAT! 8F10.5) 

11  =  1 

IF(IFLAGl.EQ.l)  GO  TO  51 

GO  TO  (4,5,6,7,8) .MATRIX 

4  CONTINUE 
DO  1  1=1,8 
J  =  J  +  1 

IF(J.GT.3)  K=K+1 
IF(J.GT.3)  J=l 
COEFOX( J,K)=INFO( I ) 
IC0UNT=IC0UNT+1 
IC=IC+1 

IF(  IC0UNT.EQ.99)  GO  TO  3 
1  CONTINUE 
GO  TO  2 

3  11=1 
IF(I1.LE.8)  11=11+1 
J  =  0 

K  =  l 

ICOUNT=0 

MATRIX=2 

5  IF( I1.GT.8)  GO  TO  2 
DO  9  1=11,8 

J  =  J  +  1 

IFU.GT.3)  K=K  +  1 
IF(J.GT.3)  J=l 
COEFCX( J, K)=INFO( I ) 
ICOUNT= ICOUNT+1 
IC=IC+1 

IF(  IC0UNT.EQ.99)  GO  TO  10 
9  CONTINUE 
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GO    TO    2 

10  11=1 
IFU1.LE.8)    11  =  11  +  1 

J  =  0 
K  =  l 

ICOUNT=0 
MATRIX=3 

6  IFdl  .GT.8)    GO    TO    2 
DO    11    1=11,8 
J=J  +  1 

IFU.GT.3)    K  =  K  +  1 
IFU.GT.3)     J=l 
C0EFMX1  J,K)  =  INFO( I) 
IC0UNT=IC0UNT+1 
IC=IC+1 
IFdC0UNT.EQ.99)     GO    TO    12 

11  CONTINUE 
GO    TO    2 

12  11=1 
IFdl.LE.8)     11  =  11  +  1 
J  =  0 
K  =  l 

ICOUNT=0 
MATRIX=4 

7  IF( II .GT.8)  GO  TO  2 
DO  13  1=11,8 
J  =  J  +  1 

IF(J.GT.ll)    K=K+1 
IF(J.GT.ll)    J=l 
CADSK J,K)=INFO( I ) 
IC0UNT=IC0UNT+1 
IC=IC+1 
IF(IC0UNT.EQ.121)     GO    TO    15 

13  CONTINUE 
GO    TO    2 

15  11=1 
IF(I1.LE.8)     11=11+1 
J  =  0 
K  =  l 

ICOUNT=0 
MATRIX=5 

8  IF! II. GT.8)  GO  TO  2 
DO  14  1=11,8 
J  =  J  +  1 

IF(J.GT.ll)    K=K  +  1 
IF(J.GT.ll)     J=l 
CFDSK  J  ,K)  =  INFO{  I  1 
ICOUNT=ICOUNT+1 
IC=IC+1 
IF( IC0UNT.EQ.121)     GO    TO    16 

14  CONTINUE 
GO    TO    2 

16  IFLAG1=1 
GO    TO    2 

50  IFLAG=1 
IF( (IFLAG.EQ.l) .AND.t IFLAG1 .EQ.l) )    GO    TO    52 

51  WRITE(6,104)     ICOJNT ,MATRI X , I C       - 
104FORMAT(»     • , • JOB    FA  I  LURE 8 , IX , « COUNT= • , I  6 , IX , 

lib, IX, « NO.PTS=' , 16) 
FLAG=l 
RETURN 

52  WRITE(6,105)     IC 

105  FORMAT (■  •  ,  '  JOB  COMPLETE1  ,  IX,  »  POI  NTS  ST0RED  =  ',I6) 
RETURN 
END 
FORTRAN 

C  SUBROUTINE  INTERP 

C 

c 
c 

SUBROUT INE  INTERP (ALFAC , FLA P, COEF  ,M) 

COMMON  COEFMX,CGEFCX,COEFOX 
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100 


101 
4 


DIMENSI 
OATA  CO 
ALFA=AL 
IF(ALFA 
IF(ALFA 
FORMAK 
lZXj'M^ 
IFiALFA 
IF(ALFA 
ALPHA=2 
IF(ALPH 
IF(ALPH 
IA=13 
WTA=0.0 
GO  TO  3 
I  A1=-AL 
IA=13-I 
ALPHA1= 
WTA=ALP 
GO  TO  3 
IA1=ALP 
IA=13+I 
ALPHA1= 
WTA=ALP 
IFtFLAP 
IF(FLAP 
IF(FLAP 
IFtFLAP 
FORMAT < 
IF(FLAP 
IB=1 
WTB=-{( 


ON    CO 

NDEG/ 

FAC*C 

.GT.l 

.GT.l 

•0'  ,  ' 

till 

.LT.- 

.LT.- 

.0*AL 

A. LI. 

A.GT. 


EF MX (3,33) ,C0EFCX(3,33) ,COE FOX ( 3 ,33 ) 

57.295773/ 

ONDEG 

0.0)     WRITE(6,100)ALFAC,M 

0.0)     ALFA=10.0 

ERROR. ALFA    EXCEEDED    L IMI TS. ALFA=' ,E 14 . 8 , 

6.0)  WRITE(6,100)     ALFACM 

6.0)  ALFA=-6.0 

FA 

0.0)  GO    TO    1 

0.0)  GO    TO    2 


PHA 

Al 

-IA1 

HA-ALPHA1 


HA 

Al 

IA1 

HA- 

.LT 

.LT 

.GT 

.GT 

«0« 

)4, 


FLAP 


ALPHA1 

.-5.0)     WRITE(6,101) 

.-5.0)     FLAP=-5.0 

.5.0)     WRITE(6f101)     FLAP 

.5.0)     FLAP=5.0 

, 'ERROR. FLAP    EXCEEDED    LIMITS. 

5,6 


FLAP=» ,E14.8) 


GO  TO 

5  IB=2 
WTB=0„0 
GO  TO  7 

6  fB=3 
WTB=<5. 

7  IF(WTA) 

8  IA2=IA- 
GO  TO  1 

9  IA2=IA 
GO  TO  1 

10  IA2=IA+ 

11  IF(WTB) 

12  IB2  =  2 
GO  TO  1 

13  IB2=IB 

14  GO  TC(1 

15  C0EF1=C 
C0EF2=C 
GO  TO  1 

16  C0EF1=C 
CUEF2=C 
GO  TO  1 

17  C0EF1=C 
C0EF2=C 

18  C0EF3=l 
GO    TO    ( 

19  C0EF1P= 
C0EF2P= 
GO  TO  2 

20  C0EF1P= 
C0EF2P= 
GO  TO  2 

21  C0EF1P= 
C0EF2P= 

22  C0EF3P= 
COEF=(C 
RETURN 
END 

FORTRAN 


5.0+FLAP)/5.0) 


O-FLAPJ/5 

8,9,10 

1 

1 


1 
1 

12, 


13,12 


5,16,17)  , 

OEFMX( IB, 

OEFMX( 

8 

OEFCX( 

OEFCX( 

8 

OEFOX( 

OEFOX( 


IB2 

IB, 
132 


M 

IA) 

,IA) 

IA) 
,IA) 


IB, 
IB2 


IA) 
,  IA) 

COEF2-C0EF1 J*ABS( WTB)+C0EF1 
19,  20,21 ) t M 
COEFtfX( I3,IA2) 
IA2) 


C0EFMX(I82 

2 

C0EFCX1 IB, 

COEFCXi IB2 

2 

COEFOX( 13, 

COEFOX( IB2 


IA2) 
,  IA2) 


IA2) 

,  I  A2  ) 
(C0EF2P-C0EF1P) *ABS ( WTB) +COEF 1 P 
0EF3P-C0EF3)*ABS(WTA)+C0EF3 
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C  SUBROUTINE    INTRP1 

C 

c 
c 

SUBROUTINE     INTRPHBETA1  f  SOB  ,COEFS  ,  M,  THE  T  A  ,  TFIETA1  ) 
COMMON    JUNK,CFDS1,CADS1 

DIMENSION    JUNK(297) ,CFDS1(11,11)  ,CADS1(  11,11) 
DATA    CONDEG/57. 295773/ 
SUB=SOB 

BETA=BETAl*CONDEG 

IF(BETA.GT.5.0)     W*ITE(6,100)    BETA1.M 
IF(BETA.GT.5.0)     BETA=5.0 
100    FORMAT (  '0«  , 'ERROR.  BETA    EXCEEDS    L I  MI TS. BETA= • , E14 . 8 , 2X, 
l»M='i  ID 
IF(BETA.LT.-5.0)     WRITE(6,100)     BETA1,M 
IF(BETA.LT.-5.0)     BETA=-5.0 
IF(BETA.LT.O.O)     GO    TO     1 
IF(BETA.GT.O.O)     GO    TO    2 
IA=6 
WTA=0.0 
GO    TO    3 

1  IA1=-BETA 
IA=6-IA1 
BETAG=-IA1 
WTA=BETA-BETAG 
GO    TO    3 

2  IA1=BETA 
IA=6+IA1 
BETAG=IA1 
WTA=BETA-BETAG 

3  IF(SUB.LT.Q.O)     SU3=0.0 
IF(SUB.GT.iO.O)     SUB=10.0 
IF(SUB.GT.O.O)     GO    TO    5 
IS=1 

WTS=0.0 
G.0  TO  7 
5  IS1=SUB 
IS=IS1+1 
SU31=IS1 
WTS=SUB-SUB1 

7  I F(WTA) 8,9,10 

8  IA2=IA-1 
GO  TO  11 

9  IA2=IA 
GO  TO  11 

10  IA2=IA+1 

11  IF(WTS.GT.O.O)  GO  TO  12 
IS2=IS 

GO  TO  15 

12  IS2=IS+1 

15  IF(MBEQ.2)  GO  TO  16 
C0EF1=CADS1(  IS, IA) 
C0EF2=CADS1( IS2,IA) 
GO  TO  18 

16  C0EF1=CFDS1( IS, IA) 
C0EF2=CFDS1(  IS2, IA) 

18  C0EF3=(C0EF2-C0EF1)*WTS+C0EF1 
IFCM.EQ.2)  GO  TO  20 

19  C0EF1P=CADS1( IS,IA2) 
C0EF2P=CADS1(IS2, IA2) 
GO  TO  22 

20  CGEF1P=CFDS1(  IS, IA2) 
C0EF2P=CF0S1( IS2, IA2) 

22  C0EF3P= (C0EF2P-C0EF1P) *WTS+CGEF1P 

C0EFS=(C0EF3P-C0EF3)*ABS( WTAJ +C0EF3 

THETA1=THETA*C0NDEG 

RETURN 

END 
//L.SYSLIB  OD       DISP=SHR,DSN  =  SYS1  .DSL  , 
//      V0L=SER=MARY,UNIT=2314 
//      DD   DISP=SHR,DSN=SYS1.F0RTLIB 
//      DD   DISP=SHR,DSN=SYS1.MPSLIB 
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//G.FT06F001  OD  SP AC£= ( C YL , ( 3 ) ) 

//G.FT03F001  DO  DSN=S1 269 .GLOP , UN IT  =  2 32 1 ,01 SP= SHR , 

//      VOL=SER=CELG04,LABEl=( ,, , IN) 

//PLOT.SYSIN  DO  * 

STEP  RESPONSE  OF  PITCH  ANGLE 
NANAGARA 
0.0        2.0        -2.0       1.         5.0        6.0 
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